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1 .  INTRODUCTION 


This  report  presents  the  results  of  studies  on  recorded  ground  acceler¬ 
ations  for  selected  rock  and  soil  sites  located  close  to  the  generating  fault 
during  two  California  earthquakes.  The  two  earthquakes  are: 

-  San  Fernando  earthquake  of  February  9,  1971  (M^  =  6.3) 

-  San  Francisco  earthquake  of  March  22,  1957  (M^  =  5.3) 

where  the  values  of  local  magnitude,  M^  are  those  given  by  Kanamori  and 
Jennings  (1978). 

The  development  of  modern  earthquake  engineering  has  included  a  signi¬ 
ficant  number  of  studies  of  recorded  strong  motion  accelerograms.  Many  of 
these  studies  have  focused  on  the  characteristics  of  the  response  spectrum 
and  on  the  value  of  the  peak  acceleration.  Studies  of  the  time  domain  char¬ 
acteristics  of  recorded  ground  motions,  such  as  their  duration  and  the  sequence 
of  the  acceleration  pulses  have  been  relatively  scarce.  However,  more  re¬ 
cently,  these  time  characteristics  are  receiving  increased  attention,  in  re¬ 
cognition  of  their  importance  to  engineering  problems  such  as  nonlinear 
response  and  liquefaction  potential  of  soil  deposits,  permanent  displacements 
of  earth  slopes,  and  the  damage  and  collapse  potential  of  yielding  structures. 
Examples  of  time  domain  studies  published  in  the  last  few  years  are  Bolt  (1973), 
Seed  et.  al .  (1975a),  Trifunac  and  Brady  (1975),  Franklin  and  Chang  (1977), 
Vanmarcke  and  Lai  (1977)  and  Dobry  et.  al .  (1977,  1978,  1978a). 

The  present  study  focuses  on  these  time  domain  characteristics,  and  uses 
recently  developed  analytical  procedures  for  processing  accelerograms  in  the 
time  domain.  A  summary  of  the  main  features  of  the  methods  used  in  this 
study  is  presented  in  Section  2  of  this  report.  The  time  domain  procedures 
mentioned  above  have  already  proven  quite  useful  for  uncovering  hidden  features 
of  recorded  accelerograms,  and  also  for  understanding  the  relation  between 
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these  features  and  relevant  seismological  parameters  (Dobry  et.  al.,  1978, 

Bond,  1979).  Most  of  the  work  done  to  date  has  focused  on  earthquake  re¬ 
cords  obtained  at  rock  sites  and  hilly  areas,  which  very  often  appear  to 
be  simpler  than  records  obtained  at  nearby  soii  sites  and  valley  areas. 

Some  results  of  this  previous  work  are  also  reproduced  and  discussed  in 
Section  2. 

A  main  objective  of  the  present  study  is  to  develop  a  better  understanding 
of  the  time  domain  characteristics  of  accelerograms  obtained  at  soil  sites, 
as  compared  with  the  characteristics  of  records  at  neighboring  rock  sites. 

A  second  objective  is  to  evaluate  the  one  dimensional  amplification  technique 
currently  being  used  to  predict  soil  motions.  Therefore,  two  criteria  were 
used  to  select  the  earthquake  records  for  analysis:  a)  several  ground  motion 
accelerograms  should  be  available  which  were  recorded  in  the  same  general 
area  during  an  earthquake,  and  at  least  one  of  the  records  must  correspond  to 
a  rock  site,  and  b)  for  some  of  the  soil  accelerograms,  there  should  be 
enough  information  available  on  local  soil  conditions  so  that  analytical  site 
response  simulations  could  be  performed.  The  records  selected  for  this  study 
fulfill  both  criteria. 

Section  3  presents  the  results  for  one  rock  and  three  soil  accelerograms 
recorded  in  the  Pasadena  area  during  the  San  Fernando  earthquake.  Section  4 
presents  similar  results  for  one  rock  and  four  soil  accelerograms  during  the 
San  Francisco  earthquake.  All  nine  accelerograms  were  processed  using  the 
time  domain  techniques  described  in  Section  2,  and  the  results  compared  to 
study  the  possible  influence  of  geologic  and  site  conditions.  In  addition, 
site  response  analyses  were  performed  at  some  of  the  soil  sites  using  as  input 
the  neighboring  rock  records.  These  analyses  used  the  widely  utilized  shear 
beam,  equivalent  linear  analysis  technique  (Program  SHAKE,  see  Schnabel  et. 
al.,  1972).  Finally,  the  simulated  soil  motions  computed  in  the  response 
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analyses  were  processed  and  compared  with  the  recorded  motions  using  the 
same  time  domain  techniques  mentioned  above.  The  results  of  these  site 
response  studies  are  also  presented  in  Sections  3  and  4. 

All  records  used  in  the  present  study  are  corrected  accelerograms 
which  were  obtained  from  California  Insitute  of  Technology,  and  which  are 
described  in  Volume  II  of  the  series.  Strong  Motion  Earthquake  Accelerograms 
(Hudson  et.  al.,  1971). 
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2.  METHODS  FOR  ANALYZING  ACCELEROGRAMS 

The  recent  interest  in  the  time  domain  characteristics  of  strong  ground 
motion  accelerograms  has  created  the  need  for  new  methods  of  analysis.  This 
need  will  no  doubt  increase  even  further  as  more  records  become  available 
from  future  earthquakes.  Therefore,  a  significant  effort  has  already  been 
spent  at  Rensselaer  Polytechnic  Institute  (R.P.I.)  in  developing  these  new 
procedures.  A  detailed  discussion  of  these  methods  and  of  their  application 
to  rock  accelerograms  is  outside  the  scope  of  this  report,  and  can  be  found 
elsewhere  (Bond,  1979).  A  summary  description  of  basic  concepts  and  procedures, 
needed  for  this  report,  follows. 

A  strong  ground  motion  accelerogram  usually  consists  of  the  three  com¬ 
ponents  of  the  particle  ground  acceleration,  a^Ct),  a^(t),  a^Ct),  where  a^ 
and  a^  are  orthogonal  horizontal  components  and  a^  is  the  vertical  component. 

The  time  t  is  the  same  for  the  three  components  and  is  measured  from  the  mo¬ 
ment  the  instrument  is  triggered.  The  triggering  time  is  not  known  on  an 
absolute  time  scale  and  it  usually  varies  between  stations.  Typical  horizontal 

components  a  (t)  and  a  (t)  for  the  San  Fernando  earthquake  are  presented  in 
x  y 

Fig.  1. 

The  procedure  used  for  the  time  domain  analysis  of  the  accelerogram  is 
to  study  several  graphs  showing  the  variation  with  time  of  the  main  character¬ 
istics  of  the  record.  These  characteristics  are:  i)  the  level  of  shaking  or 
energy  buildup  of  the  motion,  ii)  the  direction  of  the  motion,  and  iii)  the 
frequency  content  of  the  motion. 

Graphs  of  these  time  domain  characteristics  are  shown  in  Figures  1  and 
3  through  7.  These  will  be  discussed  in  more  detail  subsequently.  The  basic 
concepts  used  in  developing  these  graphs  are  those  proposed  originally  by 
Arias  (1969),  Husid  (1973)  and  Penzien  and  Watabe  (1975)  and  later  extended 
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Figure  1  -  HORIZONTAL  (ax  and  ay)  ACCELERATION  COMPONENTS,  AND 
VARIATION  OF  RMS  FOR’ THE  VECTOR  ACCELERATION,  a^ 
WITH  TIME,  445  FIGUEROA  STREET 
(Bond,  1979) 


by  Dobry  et.  al .  (1978a)  and  Bond  (1979). 


where  =  total  duration  of  accelerogram.  Therefore,  a  normalized  Husid 
plot  goes  from  h  =  0  at  the  beginning  to  h  =  1  (or  100%)  at  the  end  of  the 
record.  Figure  3  presents  the  normalized  Husid  plot  for  the  acceleration 
components  shown  in  Fig.  1,  and  recorded  at  445  Figueroa  Street  (Station 
C054  in  the  map  of  Fig.  2).  Figure  6  compares  the  normalized  Husid  plots 
for  ^igueroa  with  those  of  four  sites  labelled  "rock  stations"  in  Fig.  2. 
(The  names  of  all  stations  included  in  Fig.  2  are  shown  in  Table  1.)  These 
five  stations  are  all  located  in  the  hilly  area  between  San  Fernando  Valley 
and  Los  Angeles  Basin,  with  the  ground  conditions  consisting  of  either 
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Terrace  deposits  and  upper  Miocene 
to  Pleistocene  sedimentary  rocks 


Rock  (Sedimentary  and  volcanic 

rock  and  basement  complex) 


£susana 


Site  Conditions  of  Stations 
#  Rock  Site 
■  Stiff  Site 


££  C054  C af.  Tech  Reference  Number 

Note:  Station  locations  are  approximate 


SAN  FERNANDO  VALLEY 
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Figure  2  -  LOCATION  MAP  OF  GROUND  RECORDING  ACCELEROGRAMS,  SAN  FERNANDO 
EARTHQUAKE,  FEB.  9,  1971 


Table  1 


Accelerograph  Stations  and  Numbers  used  in  Map  of  Fig.  1 
San  Fernando  Earthquake 


CIT  No.  (*) 

C054 

E078 

1137 

K157 

0198 

L166 

G106 

G107 

G108 

G110 


Station 

445  Figueroa 

Water  and  Power  Bldg. 

15910  Ventura  Blvd. 

420  South  Grand  Ave. 

Griffith  Observatory 
3838  Lankershim 
Seismological  Lab. 

Athenaeum  (Caltech  Campus) 
Millikan  Library  (Caltech  Campus) 
Jet  Propulsion  Laboratory 


(*)  CIT  =  California  Institute  of  Technology 
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Figure  3  -  NORMALIZED  HUSID  PLOT  AND  VARIATION  OF  RMS  WITH  TIME  FOR 
RECORDED  GROUND  ACCELERATIONS,  445  FIGUEROA  STREET 
(Bond,  1979) 
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sure  4  -  VARIATION  OF  RMS  AND  FOR  RECORDED  GROUND  ACCELERATIONS, 
445  FIGUEROA  STREET  (Bond,  1979) 
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Figure  6  -  NORMALIZED  HUSID  PLOTS  FOR  RECORDS  AT  "ROCK"  SITES,  SAN  FERNANDO  EARTHQUAKE 
(Bond,  1979) 
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Figure  7  -  NORMALIZED  HUSID  PLOTS  AND  VARIATION  OF  RMS  WITH  TIME  FOR  ACCELERATIONS 
RECORDED  AT  "ROCK"  AND  SOIL  SITES  (time  scale,  t  and  t  values  are 
from  Lankershim)  (Bond,  1979)  1  2 


rock  or  shallow  deposits  (less  than  200  ft.)  of  stiff  soil  overlying  rock 
(Grant  et.  al.,  1978).  Figure  7  compares  the  normalized  Husid  plots  for  re¬ 
cords  obtained  at  a  "rock"  site  (Lankershim,  Station  L166  in  Fig.  2)  and  at  a 
soil  site  in  the  San  Fernando  Valley  (Ventura,  Station  1137  in  Fig.  2) . 

It  has  been  suggested  (Dobry  et.  al.,  1978a)  that  the  shapes  of  Husid 

\ 

plots  for  records  obtained  at  many  rock  stations  during  small  and  moderate 
shallow  earthquakes  indicate  the  existence  of  three  parts  in  the  record. 

These  are  correlated  with  three  different  sets  of  wave  arrivals  to  the  site, 
as  follows: 

1)  Initial  weak  part,  corresponding  to  P-wave  arrivals,  and  to  the  initial 
flat  slope  of  the  Husid  plot. 

2)  Intermediate  strong  part  between  times  t^  and  t^,  corresponding  mainly 

to  the  arrival  of  S-waves  following  a  direct  path  between  earthquake 
source  and  recording  station.  This  strong  part  corresponds  to  the  line 
with  a  steep  slope  of  the  Husid  plot,  it  has  a  duration  A  =  ~  t^  and 

can  be  approximated  for  engineering  purposes  by  a  straight  line. 

3)  Final  weak  part  after  time  t^,  corresponding  to  the  arrival  of  indirect 
body  waves  and  surface  waves. 

The  hypothesis  that  the  time  t^,  marking  the  beginning  of  the  steep  part  of 
the  Husid  plot  corresponds  to  the  first  S-wave  arrival,  can  be  verified  by 
independent  evidence  for  the  San  Fernando  records  shown  in  Figs.  3  through  7. 
For  these  and  other  San  Fernando  accelerograms,  Berrill  (1975)  determined  the 
first  S-wave  arrival  using  the  technique  described  by  Hanks  (1975),  which  re¬ 
quires  the  computation  of  a  set  of  rotated  horizontal  accelerograms  parallel 
and  perpendicular,  respectively,  to  the  source-station  direction.  For  the 
Figueroa  station,  Berrill  gives  t^  =  1.8  seconds,  and  this  value  has  been 
superimposed  on  Figure  3.  A  aamparison  between  the  values  of  t^  based  on 
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the  shape  of  the  Husid  plot,  and  those  given  independently  by  Berrill,  is 
presented  in  Table  2  for  several  records.  Good  agreement  exists  between 
the  two  methods.  All  of  this  confirms  that  the  beginning  of  the  strong, 
steep  part  in  the  Husid  plot  can  be  used  for  a  visual  estimate  of  the  time 
of  first  S-wave  arrival,  for  both  rock  and  soil  records. 

The  duration  of  direct  S-wave  arrivals,  (strong  part)  has  been  esti¬ 
mated  to  be  A  -  6  seconds  for  the  San  Fernando  earthquake  records  considered, 
and  the  corresponding  end  time  of  the  strong  part,  t^  =  t^  +  6  has  also  been 
included  in  Figs.  1,  3  and  6.  It  should  be  noticed  that  t^  generally  coin¬ 
cides  in  these  figures  with  the  time  at  which  the  Husid  curves  become  flatter. 
Furthermore,  the  similarity  between  all  curves  in  Fig.  6  shows  that  the  energy 
buildup  with  time  was  quite  consistent  for  ground  motions  over  the  hilly  area 
South  and  SE  of  the  San  Fernando  Valley. 

In  addition  to  the  Husid  plot,  a  second  graph  is  also  used  to  study  the 
variation  of  the  level  of  shaking  with  time.  This  is  the  RMS(t)  (Root-Mean- 
Square)  horizontal  vector  acceleration,  computed  for  a  moving  time  window  of 
width  6  (seconds).  RMS(t)  is  defined  as: 

t+0.56 


RMS (t) 


(4) 


RMS(t)  is  also  independent  of  the  orientation  of  the  horizontal  components  x 
and  y.  The  value  of  RMS(t)  is  closely  related  to  the  slope  of  the  Husid  plot 
at  the  same  time  t.  Therefore,  it  should  be  expected  that  the  values  of  RMS(t) 
are  larger  during  the  strong  part  of  the  record  than  for  t  <  t^  or  t  > 

This  is  illustrated  by  the  plot  of  RMS(t)  versus  t  presented  in  Fig.  3  for  the 
Figueroa  records  computed  using  6  =  0.5  seconds.  Figure  3  shows  the  corres¬ 
pondence  between  the  RMS(t)  and  the  Husid  plot  for  the  same  station.  The  RMS(t) 
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graph  has  proven  extremely  useful  to  study  the  sequence  of  the  main  acceler¬ 
ation  pulses,  and  also  to  understand  them  in  terms  of  wave  arrivals  to  the 
station. 

Dobry  et.  al .  (1979)  and  Bond  (1979)  have  demonstrated  that  the  particular 
shape  of  RMS(t)  for  t  <t<t2  illustrated  by  Fig.  3  is  essentially  the  same  for 
rock  and  soil  stations  covering  a  large  area  South  and  SE  of  the  generating 
fault  during  the  San  Fernando  earthquake;  these  include  all  stations  plotted 
in  Figs.  6  and  7.  What  this  says  is  that  the  peaks  of  the  RMS(t)  graph  do 
not  occur  at  random  but  correspond  to  definite  wave  arrivals,  which  are  the 
same  for  different  stations.  These  wave  arrivals  should  be  correlated  with 
the  source  mechanism  and  the  propagation  characteristics  of  the  seismic  waves 
between  the  source  and  the  station.  Furthermore,  the  horizontal  peak  acceler¬ 
ations  of  the  records  occur  much  more  often  at  the  peaks  of  the  RMS(t)  graph 
than  at  the  valleys  of  the  same  graph.  This  is  illustrated  by  Fig.  7,  which 
compares  the  Husid  and  RMS(t)  graphs  between  a  "rock"  and  soil  station 
(Stations  L166  and  1137  in  Fig.  1).  The  normalized  Husid  plots  in  Fig.  7 
are  very  different,  and  indicate  a  significant  longer  duration  of  motion  at 
the  soil  site.  The  reason  for  this  becomes  clear  by  a  comparison  of  the 
RMS(t)  graphs  in  the  same  figure.  For  t^  <  t  <  t^,  the  two  RMS(t)  curves  are 

very  similar,  and  have  the  same  sequence  and  locations  of  peaks  and  valleys, 

(*) 

indicating  the  same  wave  arrivals  at  both  stations  .  However,  for  t  >  t^, 
the  "rock"  Lankershim  record  essentially  ends,  with  RMS ( t )  dropping  below 
0.05  g.  On  the  other  hand,  for  t  >  t^  the  soil  record  has  RMS(t)  values 
larger  than  0.05  g  for  an  additional  10  seconds.  These  motions  after  the 

(*)  The  two  main  wave  arrivals  occur  at  about  t^  +  3  and  t;  +  5  as  indicated 
by  the  major  peaks  in  RMS(t)  in  both  Figs.  3  and  7.  The  peak  accelera¬ 
tions  of  the  horizontal  components  and  of  the  vector  accelerations,  axv 
also  occur  at  about  these  times  (also  see  Dobry  et.  al.,  1979). 
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end  of  the  strong  part  in  the  soil  record  can  be  explained  by  surface  waves 
trapped  in  the  alluvium  of  the  San  Fernando  Velley,  which  did  not  propagate 
into  the  rock  of  the  hills  south  of  the  valley  due  to  the  high  impedance 
ratio  between  rock  and  alluvium. 

Throughout  the  rest  of  this  report,  normalized  Husid  plots  and  RMS(t) 
graphs  are  shown  for  all  nine  records  considered.  For  the  RMS(t)  calculations, 
a  value  5  =  0.5  sec.  was  generally  used. 


2 . 2  Direction  of  Ground  Motion 

It  is  possible  to  define  a  "ground  acceleration  tensor",  [G ( t ) ]  at  a 
given  station  for  a  particular  earthquake: 


[  G  ( t )  ]  = 


XX 

8xy 

8XZ 

yx 

8yy 

8yz 

'zx 

00 

8zz 

t+0.56 


where  g . ,  =  g. . (t)  = 
ij  ij 


f  a. (t)  a. 

J  1  J 


t-0.56 

of  the  matrix  have  the  form: 

t+0.56 


g  (t)  = 
Bxx 


•/ 


a  (t)dt 
x 


t-0.56 

while  the  off-diagonal  terms  have  the  form: 

t+0.56 


g  (t)  =  g 
sxy  vx 


/ 


(5) 


.  (t)dt.  Therefore,  the  diagonal  terms 


a  (t)  a  (t)dt 
x  y 


t-0.55 
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As  before,  x,  y  are  horizontal  axes  and  z  denotes  the  vertical  axis. 

The  symmetric  tensor  [ G ( t ) ]  varies  with  time,  and  it  is  also  a  function 
of  the  window  width,  6.  Arias  (1969)  and  Penzien  and  Watabe  (1975)  have  shown 
that  [G ( t ) ]  is  indeed  a  tensor,  as  its  components  for  any  new  set  of  orthogonal 
axes  x' ,  y' ,  z'  can  be  computed  by  using  the  rules  of  tensor  transformation. 

It  is  possible  to  compute  the  principal  directions  of  this  acceleration 
tensor  using  standard  procedures.  Of  special  interest  is  the  angle,  ^  between 
the  major  principal  direction  (major  principal  plane)  and  the  vertical  axis. 

Some  results  by  Kubo  and  Penzien  (1976)  using  a  very  large  value  of  6  =  5  seconds 
suggested  that  is  close  to  0  for  t  <  t^  and  that  it  switches  to  =  90°  for 
t  >  t^.  This  means  that  the  acceleration  is  predominantly  vertical  (P-wave)  for 
t  <  t^  and  predominantly  horizontal  (SH-wave)  for  t  >  t^.  That  this  is  indeed 
the  case,  and  the  fact  that  the  <f^(t)  graph  calculated  with  a  much  narrower  win- 
down,  6  =  0.5  sec.  can  be  very  helpful  in  determining  t^  and  for  analyzing 
ground  accelerograms,  is  shown  conclusively  by  Bond  (1979).  This  is  illustrated 
by  the  plot  of  ji,  versus  t  for  the  Figueroa  record  in  Fig.  4.  Figure  4  compares 

•L 

the  <f^(t)  and  RMS(t)  graphs  for  the  same  record,  and  clearly  demonstrates  that 
these  two  plots  supplement  each  other.  Of  particular  interest  is  the  jump  in 
both  graphs  at  t^  -  1.8  seconds. 

In  this  report,  j>^(t)  graphs  were  systematically  plotted  for  6  =  0.5 
seconds  and  are  discussed  throughout  the  text.  The  corresponding  $^(t)  graphs, 
which  give  the  angle  between  the  minor  principal  direction  (minor  principal 
plane)  and  the  vertical  axis,  were  also  obtained  and  are  included  in  Appendix  1. 
The  angle  $  (t),  and  the  azimuth  angles  with  respect  to  the  North  of  the  three 
principal  directions,  9^(t),  and  8^(0  can  also  be  computed  and  plotted 

versus  t.  In  general,  is  the  most  important  of  all  these  angles.  The  use 
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of  graphs  for  some  of  the  other  angles  in  the  analysis  and  interpretation 
of  accelerograms  is  discussed  by  Bond  (1979). 


2.3  Frequency  Content  of  the  Motion 


The  graphs  discussed  in  Sections  2.1  and  2.2  are  very  useful  for  time- 
domain  analysis  of  accelerograms.  They  depict  in  a  simple  way  the  variation 
with  time  of  level  of  shaking  and  energy  buildup  of  the  record,  and  of  the 
direction  of  the  ground  accelerations.  Furthermore,  these  graphs  are  indepen¬ 
dent  of  the  orientation  of  the  horizontal  components  of  the  instrument,  which 
eliminates  the  need  to  consider  this  additional  variable.  It  would  be  desirable 
to  have  a  similar  simple  graph  for  the  predominant  frequency  of  the  ground 
acceleration,  f(t)  (cps),  which  could  provide  information  on  changes  in  the 
frequency  content  associated  with  different  parts  of  the  record  or  with 
specific  wave  arrivals.  It  would  also  be  convenient  that  this  graph  be 
independent  of  the  orientation  of  the  instrument. 

There  are  several  definitions  of  predominant  frequency  of  the  motion 
(i.e.,  Crandall  and  Mark,  1963;  Clough  and  Penzien,  1975).  Also,  different 
techniques  have  been  developed  for  characterizing  and  plotting  the  time  varia¬ 
tion  of  the  frequency  content  of  the  motion,  which  are  mostly  based  on  the 
use  of  multi-filtering  procedures  (Liu,  1970  and  1971;  Trifunac,  1971;  Perez, 
1973;  Saragoni  and  Hart,  1973;  Kameda,  1975;  Joannon  et.  al.,  1977).  These 
and  other  definitions  and  techniques  are  presently  being  studied  at  RPI,  with 
the  purpose  of  developing  the  desired  graph  of  predominant  frequency  versus 
time . 


For  the  purpose  of  this  report,  a  crude  technique  was  used  to  estimate 


f(t).  This  consisted  of  counting  the  number  of  peaks,  N,  of  the  vector 

2  2  2 

horizontal  acceleration  squared,  a  (t)  =  a  (t)  +  a  (t),  between  the  times 

xy  x  y 

t  -  0.5tSj  and  t  +  0.56,,  and  defining  f(t)  by  Eq.  6: 
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(6) 


f(t) 


N 


The  graph  of  f(t)  versus  t  for  the  Figueroa  record  is  presented  in  Fig.  5, 
for  the  time  window  <5^  =  1  sec.  Figure  5  also  reproduces  the  corresponding 
graph  of  RMS(t).  The  values  of  f(t)  in  Fig.  5  during  the  strong  part  are 
between  4  and  6  cps.  A  similar  value  for  the  average  f(t)  in  the  range 
t^  <  t  <  t^  was  found  for  San  Fernando  records  by  Dobry  et.  al.  (1979)  using 
a  more  refined  computational  method. 

Graphs  of  f(t)  using  Eq.  (6)  and  =  1  second  were  obtained  for  all 
records  considered  in  this  study,  and  are  included  in  Appendix  2. 
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3.  GROUND  MOTIONS  IN  THE  PASADENA  AREA 


(San  Fernando  Earthquake,  February  9,  1971) 

The  San  Fernando,  California  earthquake  of  February  9,  1971  was  a  thrust 
fault  event  occurring  along  the  San  Fernando  Fault.  The  fault  ruptured  along 
an  estimated  15-23  km  length,  and  it  broke  to  the  ground  surface  at  the  North 
end  of  the  San  Fernando  Valley  as  shown  in  Fig.  2.  Estimates  for  the  total 
duration  of  the  fault  rupture  vary  between  6  and  9  seconds  (Dobry  et.  al., 
1978a).  Values  of  local  magnitude,  M^  =  6.3  and  surface  wave  magnitude. 

Mg  =  6.6  have  been  reported  by  Kanamori  and  Jennings  (1978). 

The  earthquake  occurred  near  Los  Angeles,  in  an  area  densely  instrumented 
for  strong-motion  measurements.  Many  useful  ground  motion  accelerograms  were 
obtained,  which  have  been  and  still  are  being  studied  by  many  workers  (i.e., 
Bolt,  1973,  Boore  and  Zoback,  1974;  Bouchon  and  Aki,  1977;  Hanks,  1975,  1976; 
Crouse,  1973;  Berrill,  1975;  Seed  et.  al.,  1974,  1975).  Time  domain  and  dura¬ 
tion  studies  performed  by  the  authors  on  San  Fernando  accelerograms  are  re¬ 
ported  by  Dobry  et.  al.  (1977,  1978a,  1979)  and  Bond  (1979).  Basic  informa¬ 
tion  on  various  aspects  of  this  earthquake  can  be  found  in  the  Report  on 
San  Fernando,  California  Earthquake  of  Feb . 9 ,  1971 ,  3  vol  ■  (Murphy,  1973). 

Four  ground  accelero-'  ^ms  recorded  in  the  Pasadena  area  were  selected  for 
this  study.  These  correspond  to  the  following  stations  in  Figs.  2  and  8. 


Fig.  2 

Fig.  8 
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Detailed  descriptions  and  discussions  of  these  accelerograms,  of  their 
response  spectra  and  of  additional  seismoscope  ground  records  from  the  same 
area  have  been  published  by  Hudson  (1972),  Jennings  (1973)  and  Berrill  (1975). 
Additional  information  on  site  conditions  at  some  of  the  stations  including 
measured  shear-wave  velocities  of  the  subsoil  are  available  in  Eguchi  et.  al . 
(1976),  and  more  recently  in  Shannon  and  Wilson  (1978). 

Figure  8  is  a  map  of  the  Pasadena  area,  originally  published  by  Gutenberg 
(1957),  and  including  the  locations  of  strong-motion  accelerographs  and 
seismoscopes  which  recorded  during  the  San  Fernando  earthquake.  The  four 
accelerograph  stations  are  SL,  ATH,  ML  and  JPL.  Additional  information  on  the 
four  stations  is  provided  in  Table  3. 

3. 1  Recorded  Motions 

The  accelerograms  at  stations  SL,  JPL,  ATH  and  ML  were  studied  using  the 
time  domain  techniques  described  in  Section  2.  As  indicated  by  Hudson  (1972), 
and  illustrated  by  the  maps  in  Figs.  2  and  8,  the  distances  between  the  sta¬ 
tions  and  the  earthquake  source  is  such  that  distance  attenuation  effects 
should  be  a  minor  factor  in  explaining  differences  between  ground  motions  at 
the  four  sites. 

Table  4  summarizes  the  values  of  the  recorded  peak  accelerations  and 
their  times  of  occurrence  for  the  NS  and  EW  components.  All  values  of  peak 
acceleration  are  in  the  range  0.1  -  0.2  g,  and  they  do  not  show  any  clear 
influence  of  site  conditions.  The  peak  accelerations  at  ML  are  about  twice 
than  those  at  ATH,  although  the  two  instruments  were  within  1,000  feet  of 
each  other  and  on  essentially  the  same  soil  conditions.  Berrill  (1975)  sug¬ 
gested  that  the  stronger  motions  at  ML  may  have  been  caused  by  standing  waves 
formed  in  the  soil  between  the  basement  of  the  ML  building  and  those  of  adja¬ 
cent  buildings.  If  that  is  the  case,  the  ATH  record  would  be  more  representative 
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of  free  field  accelerations  than  the  ML  record. 


Another  interesting  fact  in  Table  4  is  the  consistency  of  the  times  of 
occurrence  of  the  peak  accelerations.  For  seven  out  of  the  eight  horizontal 
components,  the  times  (t^  -  t^)  are  clustered  around  the  narrow  range  2.6  - 
3.4  seconds,  and  for  one  component,  (t^  -  t^)  =  2.2  seconds.  This  is  con¬ 
sistent  with  the  discussion  of  Figs.  3  and  7  in  Section  2,  and  suggests  that 
the  wave  arrival  at  t  ~  t^  +  3  seconds  was  more  prominent  in  the  Pasadena 
area  than  the  second  arrival  at  t  «  t,  +  5  sec. 

The  results  of  the  time  domain  analyses  of  the  accelerograms  recorded  at 
the  four  stations  are  presented  in  Figs.  9  through  18.  These  include  normalized 
Husid  plots,  and  RMS(t)  and  $^(t)  graphs.  (The  <f>j(t)  and  f(t)  graphs  are  in¬ 
cluded  in  Appendices  1  and  2.)  In  al'  these  figures,  the  estimated  values  of 
t^  from  Berrill  (1973)  and  of  t^  =  t^  +  6  have  also  been  superposed.  Figure 
19  includes  the  acceleration  response  spectra  for  the  SL  (rock)  and  ATH  (soil) 
stations . 

Figures  9  and  10  include  the  normalized  Husid  plots  for  the  four  stations. 
For  all  four  plots,  the  locations  of  t^  given  by  Berrill  are  consistent  with 
the  shape  of  the  curves.  The  comparison  in  Fig.  9  indicates  a  rate  of  energy 
buildup  almost  identical  for  SL  and  JPL,  with  h(t^)  -  83  to  90%.  Figure  10 
indicates  a  rate  of  energy  buildup  very  similar  for  the  two  nearby  ATH  and  ML 
soil  stations,  with  h ( t ^ )  “  70  to  75%.  The  comparison  with  the  SL  plot  in  Fig. 
10  also  suggests  that  some  additional  wave  reflections  occurred  after  t^  in 
the  soil  area  represented  by  ATH  and  ML  which  were  not  present  on  rock  at  SL. 
This  is  similar  (although  not  as  prominent)  to  the  wave  reflections  in  the  San 
Fernando  Valley  soil  stations  illustrated  by  Fig.  7.  It  also  becomes  apparent 
if  the  Trifunac  and  Brady  (1975)  duration  D  for  the  four  records  are  compared. 

(D  is  defined  as  the  time  needed  to  buildup  the  ordinates  in  the  Husid  plot 
from  5%  to  95%. )  The  values  of  D  are: 
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=  5.0  sec. 


Station 


Duration,  D 


SL  (rock) 

9. 

.  1  sec 

JPL 

10. 

,6 

ATH 

13 

ML 

15. 

,6 

These  longer  durations,  D  for  motions  at  some  soil  sites  as  compared  with 
rock  are  consistent  with  similar  results  for  the  San  Fernando  earthquake  and 
for  other  earthquakes  found  by  Trifunac  and  Brady  (1975)  and  Dobry  et.  al. 

(1977,  1978a). 

In  Fig.  11  the  Husid  plots  at  SL  and  at  one  of  the  "rock"  sites  included 
in  Fig.  6,  are  plotted  together.  This  comparison  reveals  that  the  rate  of 
energy  buildup  at  rock  sites  in  the  Pasadena  area  and  in  other  areas  during 
the  San  Fernando  quake  was  generally  consistent. 

Figures  12,  13  and  14  present  the  RMS(t)  and  <ji^(t)  graphs  for  the  SL  (rock) 
and  JPL  (soil)  stations.  The  comparison  in  Fig.  14  indicates  that  the  RMS(t) 
graphs  are  very  similar  for  both  stations,  with  a  sharp  jump  at  t  =  t^  and  a 
larger  peak  at  t  -  t^  +  3,  as  expected.  The  shape  of  the  RMS  graphs  in  the 
strong  part  is  indeed  very  similar  between  the  two  stations,  including  both 
the  shapes  and  the  ordinate  values.  Figure  14  also  shows  a  common  wave  arrival 
at  t  =  t^  +  30  sec.  for  both  stations.  Figures  12  and  13  show  generally 
60°  <  ^  <  90°  during  the  strong  part,  with  the  shift  of  |  at  t  =  t^  from  low 

values  to  tj> ^  >  60°  being  particularly  clear  for  JPL  in  Fig.  13. 

Figures  15  through  17  present  the  RMS(t)  and  <£^(t)  graphs  for  the  two 
nearby  soil  stations  ATH  and  ML.  For  both  stations,  the  values  of  t^  coincide 
well  with  the  sharp  jumps  of  RMS  and  ^ ,  with  shifting  from  #  f  20°  at 
t  <  t^  to  60°  <  <t>^  <  90°  in  the  strong  part,  t^  <  t  <  t^.  Figure  16  reveals 
again  that  the  RMS  values  have  a  peak  at  t  =  t^  +  3,  although  the  peak  is 
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RMS  ACCELERATION(mm/sec  )  x  10  MAX.  PRINCIPAL  PLANE-PHI 


Figure  13  -  VARIATION  OF  RMS  AND  0  i  WITH  TIME  FOR  RECORDED  ACCELERATIONS, 
JET  PROPULSION  LABORATORY,  PASADENA 
(0.5  -  SEC.  WINDOW; 


t,  =1.8  sec. 


RMS  ACCELERATION(mm/sec2)  x  lO1  M AX .  PRI NC I  PAL  PL  AN E  -  P H I 


Figure  15  -  VARIATION  OF  RMS  AND  Oj  WITH  TIME  FOR  RECORDED  ACCELERATIONS, 
MILLIKAN  LIBRARY,  PASADENA 
(0.5  -  SEC.  WINDOW) 
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t2  =11.0  sec. 


Figure  16  -  COMPARISON  OF  RMS  CURVES  FOR  RECORDED  HORIZONTAL 
ACCELERATIONS,  CALTECH  CAMPUS,  PASADENA 
(0.5  -  SEC.  WINDOW) 

(time  scale,  t^  and  values  are  from  Athenaeum) 
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Figure  17  -  VARIATION  OF  WITH  TIME  FOR  RECORDED  ACCELERATIONS,  ATHENAEUM,  PASADENA 

(0.5  -  SEC.  WINDOW) 


significantly  larger  for  the  ML  curve,  as  could  be  expected  from  the  larger 
values  of  the  peak  accelerations  at  ML  (see  Table  4).  The  shapes  of  the  RMS 
curves  in  Fig.  16  are  similar,  with  the  peaks  and  valleys  for  the  two  curves 
being  generally  at  the  same  locations. 

Figure  18  compares  the  RMS (t )  graph  at  ATH  (soil)  with  that  at  SL  (rock). 
The  general  ordinates  and  the  shapes  of  both  curves  are  similar,  with  the 
major  peak  occurring  at  t  =  t^  +  3.  However,  the  agreement  in  the  details  of 
the  shapes  is  less  than  what  was  found  between  SL  and  JPL  in  Fig.  14,  and 
between  ATH  and  ML  in  Fig.  16.  In  addition,  in  Fig.  18,  the  RMS  values  for 
ATH,  which  were  somewhat  smaller  during  most  of  the  strong  part,  t  <  t^,  be¬ 
come  larger  than  those  for  SL  after  t This  is  consistent  with  the  dif¬ 
ferences  in  the  shapes  of  the  Husid  plots,  discussed  above,  and  with  the 
hypothesis  of  additional  wave  reflections  in  the  area  in  which  the  soil 
stations  ATH  and  ML  are  located. 

Figure  19  compares  the  acceleration  response  spectra  for  the  records  at 
SL  and  ATH,  for  a  5%  damping  ratio.  At  SL,  the  accelerations  and  spectral 
values  are  considerably  stronger  for  the  EW  than  in  the  NS  component,  while 
at  ATH  the  accelerations  and  spectral  values  for  NS  and  EW  are  similar.  This 
directional  aspect  of  the  SI.  records  has  been  noticed  by  Hudson  (1972). 

3 . 2  Site  Response  Simulat ions 

Site  response  calculations  were  performed  at  the  Athenaeum  site,  using 
as  rock  input  the  horizontal  motions  recorded  at  the  Seismological  Laboratory. 
The  soil  profile  used  for  the  calculations,  including  the  maximum  values  of 
the  shear  wave  velocity  at  small  strains,  Vg,  is  shown  in  Fig.  20.  The  Vg 
values  between  ground  surface  and  a  depth  of  about  25  feet  was  obtained  from 
refraction  surveys  at  the  site,  reported  by  Eguchi  et.  al.  (1976).  Below 
25  feet,  the  Vg  values  were  obtained  from  the  correlations  for  "Old  Alluvium" 
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RMS  ACCELERATION  (mm/sec  )x  10 


t,  =  5  sec. 


TIM  E  (sec.) 


Figure  18  -  COMPARISON  OF  RMS  CURVES  FOR  HORIZONTAL  ACCELERATIONS 
RECORDED  AT  ROCK  AND  SOIL  SITES,  PASADENA  AREA 
(0.5  -  SEC.  WINDOW) 

(time  scale,  and  values  are  from  Athenaeum) 
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Figure  19  -  COMPARISON  OF  RESPO 
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Figure  20  -  SOIL  PROFILE  AT  ATHENEUM  SITE  USED  FOR  SITE 
RESPONSE  SIMULATIONS. 


published  by  Campbell  and  Duke  (1976)  for  the  Los  Angeles  Area. 

Shear  wave  velocity  refraction  surveys  performed  at  the  Millikan  Library 
site  together  with  other  information  published  by  Eguchi  et.  al.  (1976),  in¬ 
dicate  that  the  ATH  and  ML  sites  are  very  similar,  and  that  the  profile  in  Fig. 

20  is  valid  also  at  ML.  Therefore,  all  simulations  performed  using  the  pro¬ 
file  in  Fig.  20  are  in  principle  applicable  to  free  field  conditions  at  both 
sites.  This  confirms  that  the  differences  between  the  recorded  motions  at 
the  two  stations,  discussed  in  Section  3.1,  were  not  caused  by  differences 
in  local  site  conditions.  Based  on  Berrill's  hypothesis  discussed  in  Section 
3.1,  it  is  believed  that  the  Athenaeum  records  are  a  better  representation 
of  the  free  field  motions.  Therefore,  the  results  of  the  simulations  were 
compared  with  the  Athenaeum  recorded  accelerograms. 

For  the  simulations,  Program  SHAKE,  developed  at  the  University  of 
California  at  Berkeley  (Schnabel  et.  al.,  1972),  was  used.  This  code  is  a 
one-dimensional  program,  which  assumes  that  all  horizontal  ground  motions  are 
shear  waves  travelling  vertically  in  the  soil.  This  seems  a  reasonable  as¬ 
sumption  for  the  strong  part,  t^  <  t  <  where  the  measured  values  of 
at  SL,  ATH  and  ML  are  close  to  90°  (see  Figs.  12,  15  and  17). 

SHAKE  assumes  that  the  soil  profile  is  linear-hysteretic  (viscoelastic)  and  it 
iterates  until  the  soil  properties  selected  (modulus  and  damping)  are  consistent 
with  the  calculated  seismic  shear  strains  in  the  soil.  For  the  simulations,  the 
standard  modulus  reduction  and  damping  curves  versus  strains  provided  in  the 
program  for  sand  soils  (see  also  Seed  and  Idriss,  1970)  were  used. 

Two  series  of  site  response  simulations  were  performed  with  SHAKE,  to 
evaluate  the  sensitivity  of  the  results  with  respect  to  the  selected  depth- 
to-rock.  In  one  of  these  series,  the  profile  shown  in  Fig.  20  was  used,  with 
a  total  soil  thickness  of  100  ft.,  and  a  value  Vg  =  2,500  fps  for  the  rock; 
this  Vg  is  consistent  with  the  near-surface  values  measured  in  the  rock  at 
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the  SL  site  (Eguchi  et.  al.,  1976).  In  the  second  series,  only  the  upper 
0-50  ft.  of  the  profile  in  Fig.  20  was  used,  with  Vg  =  2,000  fps  for  the  rock. 
In  all  calculations,  the  SL  input  motions  were  assumed  to  act  at  a  rock  out¬ 
crop. 

Table  5  presents  a  comparison  between  computed  and  simulated  peak  hori¬ 
zontal  accelerations  and  their  times  of  occurrence,  (t  -  t^).  The  comparison 
between  acceleration  values  is  reasonably  close  for  the  NS  component,  but  the 
simulated  values  are  two  to  three  times  larger  than  the  recorded  accelerations 
for  the  EW  component.  This  is  related  to  the  uni-directional  character  of  the 
input  SL  accelerograms  already  mentioned  in  Section  3.1.  On  the  other  hand, 

the  simulated  values  of  (t  -  t,)  are  closer  to  the  recorded  value  in  the  EW 

P  1 

direction. 

Figures  21  and  22  present  comparisons  between  the  simulated  and  recorded 
acceleration  response  spectra  for  the  NS  and  EW  components,  respectively.  In 
the  NS  direction  (Fig.  21)  the  simulated  spectral  values  are  larger  than  the 
recorded  values  for  periods  up  to  about  0.30  seconds,  while  the  opposite  is 
true  for  larger  periods.  This  is  similar  to  the  comparison  between  NS  re¬ 
corded  spectra  for  ATH  and  SL  in  the  upper  part  of  Fig.  19.  In  the  EW  direc¬ 
tion  (Fig.  22),  the  simulated  spectra  is  much  larger  than  the  recorded  spec¬ 
tra  for  periods  up  to  about  0.3  seconds,  while  they  are  similar  for  periods 
larger  than  0.3  seconds;  this  is  in  turn  similar  to  the  comparison  between 
recorded  ATH  and  SL  spectra  in  the  EW  direction  shown  in  the  lower  part  of 
Fig.  19.  In  summary,  the  characteristics  of  the  simulated  soil  spectra  seem 
to  be  determined  mainly  by  the  characteristics  of  the  corresponding  input 
rock  spectra. 

Time  domain  analyses  of  the  simulated  soil  accelerograms  at  the  ATH  Site 
were  also  performed.  Normalized  Husid  plots  and  RMS(t)  graphs  were  developed 
and  compared  with  the  corresponding  recorded  curves  at  the  site.  $^(t)  graphs 


Simulated  and  Recorded  Peak  Horizontal  Accelerations 


Absolute  Accl.in  G's 


- Simulated  (H  =  100  ft.) 


Period  in  sec. 


Figure  21  -  COMPARISON  OF  RESPONSE  SPECTRA  FOR  RECORDED  AND 

SIMULATED  ACCELERATIONS,  ATHENAEUM  SITE,  NS  COMPONENT 
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could  not  be  obtained  as  no  simulation  was  performed  of  the  vertical  motions 
at  the  ATH  Site. 

The  normalized  Husid  plot  for  the  simulation  using  V<.  =  2,500  fps  for  the 
rock  has  been  superimposed  on  Fig.  10.  Figure  10  shows  that  the  rate  of  energy 
buildup  of  the  simulated  soil  motions  is  identical  to  that  of  the  input  (SL) 
record.  Therefore,  the  increase  in  duration  at  the  ATH  and  ML  Sites,  discussed 
in  Section  3.1,  is  not  predicted  by  the  site  response  analyses  using  SHAKE. 

This  is  confirmed  by  a  comparison  of  the  5  to  95%  durations  D  of  the  motion 
obtained  from  the  Husid  plots  in  Fig.  10: 


Description 
SL  Crock),  recorded 
ATH  (soil),  recorded 
ATH  ,  simulated 

H  =  100  ft. 
ATH  ,  simulated 

H  =  50  ft. 


Duration,  D 
9.1  sec. 
13 
9.0 

9.0 


The  comparisons  for  the  RMS(t)  graphs  are  shown  in  Figs.  23  and  24. 

Figure  23  compares  the  curve  of  the  simulated  ground  motions  at  the  Athenaeum 
site,  with  the  corresponding  RMS(t)  for  the  rock  input.  The  curve  for  the 
simulation  has  larger  ordinate  values,  but  the  shapes  of  the  two  curves  are 
very  similar,  with  essentially  a  total  coincidence  in  the  locations  of  peaks 
and  valleys.  The  comparison  between  simulated  and  recorded  RMS  curves  at 
Athenaeum  is  shown  in  Fig.  24.  The  shape  predicted  by  the  simulations  in  the 
strong  part  is  generally  similar  to  that  of  the  recorded  curve;  in  particular, 
the  location  of  the  large  recorded  peak  at  t  5  t^  +  3  is  predicted  by  the 
simulations.  However,  the  details  of  the  shape  do  not  show  such  good 
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Figure  23  -  COMPARISON  OF  RMS  CURVES  FOR  INPUT  AND  OUTPUT  ACCELERATIONS,  SOIL  RESPONSE  SIMULATION, 
ATHENAEUM  SITE  (time  scale,  t.  and  t_  values  are  from  Seismological  Lab.) 


«  5.0  sec 
=  11.00  sec 


correspondence,  and  the  agreement  in  the  locations  of  peaks  and  valleys  is 


not  better 
18.  Also, 
corded  RMS 


than  that  between  recorded  curves  at  ATH  and  SL  stations  in  Fig. 
the  simulated  RMS  values  are  significantly  larger  than  the  re¬ 
ordinates  in  Fig.  24. 
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4.  GROUND  MOTIONS  IN  THE  SAN  FRANCISCO  BAY  AREA 
(San  Francisco  Earthquake,  March  22,  1957) 

The  San  Francisco,  California  earthquake  of  March  22,  1957  was  a  small 
event,  with  a  local  magnitude,  M^  =  5.3  (Kanamori  and  Jennings,  1978).  How¬ 
ever,  it  occurred  along  the  segment  of  the  San  Andreas  Fault  closer  to  down¬ 
town  San  Francisco,  and  some  damage  to  buildings  did  occur  (Steinbrugge  et. 
al.,  1959).  For  an  earthquake  of  this  magnitude,  it  can  be  estimated  that  the 
length  of  the  fault  rupture  did  not  exceed  a  few  kilometers  and  the  total 
duration  of  the  rupture  may  have  been  about  2  seconds  (Dobry  et.  al.,  1978a). 

Ground  motions  during  the  earthquake  were  recorded  by  strong  motion 
accelerographs  at  five  stations  located  in  San  Francisco  and  Oakland.  Studies 
of  this  earthquake  and  of  the  recorded  ground  accelerograms  include  Hudson 
and  Housner  (1958),  Oakeshott  (1959),  Idriss  and  Seed  (1968),  Dobry  et.  al. 
(1978a)  and  Singh  et.  al.  (1979).  These  studies  have  included  computation  of 
response  spectrum,  instrumental  intensity  and  duration,  as  well  as  analyses 
and  discussions  of  the  influence  of  local  site  conditions.  The  papers  by 
Idriss  and  Seed  and  Singh  et.  al.  include  information  on  site  conditions  and 
soil  properties  at  the  stations. 

Figure  25  is  a  map  showing  the  five  stations  and  the  estimated  rupture 
zone  along  the  San  Andreas  Fault.  One  of  the  instruments  was  located  on  rock 
in  Golden  Gate  Park,  while  three  of  the  stations  were  soil  sites  in  downtown 
San  Francisco.  The  soil  profile  for  these  downtown  sites  consists  of  sand  and 
clay  deposits  ranging  in  thickness  from  140  to  285  feet.  One  of  these 
(Southern  Pacific  Building),  is  a  clay  site,  and  has  been  classified  as  a 
"soft  to  medium"  site  from  the  viewpoint  of  the  influence  of  local  conditions 
on  earthquake  ground  motions  (Seed  et.  al . ,  1974,  1975).  The  fifth  instrument 
was  located  across  the  San  Francisco  Bay,  in  Oakland,  on  1,000  feet  of  soil, 


Figure  25  -  LOCATION  MAP  OF  GROUND  RECORDING  ACCELEROGRAPHS , 
SAN  FRANCISCO  EARTHQUAKE,  MARCH  22,  1957 
(Idriss  and  Seed,  1968) 
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with  the  upper  100  feet  consisting  mainly  of  dense  sands.  Some  additional 
information  on  the  five  stations  is  provided  in  Table  6. 

4 . 1  Recorded  Motions 

Table  7  summarizes  the  values  of  the  recorded  peak  accelerations  and  of 
their  times  of  occurrence  for  all  horizontal  components  at  the  five  stations. 
The  peak  accelerations  range  from  about  0.10  g  at  Golden  Gate  Park  to  0.02- 
0.04  g  at  Oakland  City  Hall.  Table  7  also  includes  the  estimated  times,  t^ 
of  the  first  S-wave  arrival. 

The  results  of  the  time  domain  analyses  of  the  accelerograms  recorded  at 
the  five  stations  are  presented  in  Figs.  26  through  32.  These  include  nor¬ 
malized  Husid  plots  as  well  as  graphs  of  RMS(t)  and  4>^(t).  The  <t>j(t)  and 
f(t)  graphs  are  included  in  Appendices  1  and  2.  Figure  33  includes  a  compari¬ 
son  between  recorded  acceleration  response  spectra  at  the  rock  site  and  at 
the  soft  Southern  Pacific  Bldg.  site. 

Figures  26  and  27  present  the  normalized  Husid  plots  for  the  five  sta¬ 
tions.  The  three  curves  for  soil  sites  in  the  downtown  area  in  Fig.  27  in¬ 
dicate  an  aftershock  at  about  30  seconds.  This  aftershock  could  also  be 
clearly  identified  in  the  original  accelerograms  recorded  at  Golden  Gate 
Park  (not  included  here),  but  was  not  present  in  the  more  distant  Oakland 
records.  The  shapes  of  the  Husid  plots  indicate  significantly  longer  dura¬ 
tion  for  the  motion  at  the  soil  sites  than  at  the  Golden  Gate  Park  rock  site. 
Of  special  interest  is  the  much  longer  duration  at  the  soft  Southern  Pacific 
Building  site.  This  longer  duration  is  due  to  long  period  motions  which 
occurred  in  the  soil  after  the  end  of  the  strong  part  of  the  record.  These 
long  period  motions  have  a  period  of  about  1  second  and  they  have  been  at¬ 
tributed  to  either  soil-structure  interaction  affects  (Hudson  and  Housner, 
1958)  or  one-dimensional  amplification  of  vertically  propagating  shear  waves 
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Characteristics  of  Accelerograph  Ground  Stations 


Characteristics  of  Recorded  Horizontal  Accelerations 


PERCENT  TOTAL  ARIAS  INTENSITY 


TIME  (sec.) 

Figure  27  -  NORMALIZED  HUSID  PLOTS  OF  RECORDED  HORIZONTAL 
ACCELERATIONS,  SOIL  SITES,  SAN  FRANCISCO 

DOWNTOWN  AREA  (time  scale,  t.  ,  t„  and  (t  ) 

12  la 

values  are  from  Southern  Pacific  Bldg.) 
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(Dobry  et.  al.,  1978a;  Singh  et.  al.,  1979). 

Figures  28  through  32  present  graphs  of  RMS(t)  and  $^(t)  for  the  five 
stations.  The  presence  of  the  aftershock  can  be  noted  in  all  these  figures  ex¬ 
cept  in  the  more  distant  Oakland  City  Hall,  Fig.  32.  It  is  of  interest  to  dis¬ 
cuss  the  aftershock  first,  as  the  instruments  were  still  running  when  the  after¬ 
shock  occurred  and  therefore  both  the  P-wave  and  S-wave  arrivals  were  recorded 
by  the  accelerographs .  In  Figs.  28  through  31,  the  arrival  of  the  aftershock  P- 
wave  is  clearly  marked  by  a  dip  in  to  values  less  than  20°,  indicating  almost 
vertical  accelerations.  The  aftershock  S-wave  arrival  is  associated  with  a  return 
of  4>^  to  80°  to  90°,  indicating  the  predominance  of  horizontal  accelerations. 

As  a  result  of  the  consistent  behavior  between  the  graphs  and  the  RMS 
and  Husid  plots,  the  arrival  times  of  the  aftershock  S-wave,  (t^)^  could  be 
well  defined  for  the  four  stations: 


Station 

(tl>a 

Golden  Gate  Park 

29.16  sec 

State  Building 

28.06 

Alexander  Building 

29.00 

Southern  Pacific  Building 

28.44 

These  values  of  (t  )  have  been  superimposed  in  Figs.  26  through  3],  and 
i.  Si 

have  been  used  to  align  the  Husid  plots  for  the  three  downtown  stations  in 

Fig.  27,  on  the  assumption  that  the  travel  paths  for  direct  waves  from  the 

fault  were  essentially  the  same  for  the  main  shock  and  the  aftershock. 

The  times  corresponding  to  the  first  S-wave  arrival  for  the  main  event, 

tj,  were  then  estimated  from  a  study  of  Figs.  26  through  32.  The  shapes  of 

the  Husid  plots,  and  of  the  RMS(t)  and  <t>^(t)  graphs  were  considered  for  this 

purpose,  together  with  the  values  of  (t.)  for  the  aftershock  determined  above. 

I  a 
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RMS  ACCELERATION  (mm/secr)x  10  MAX  PRINCIPAL  PLANE-PHI 


TIME  (sec) 

Figure  28  -  VARIATION  OF  RMS  AND  <p1  WITH  TIME  FOR  ACCELERATIONS  RECORDED 
AT  ROCK  SITE,  GOLDEN  GATE  PARK,  SAN  FRANCISCO 
(0.5  -  SEC.  WINDOW) 
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RMS  ACCELERATION  (mm/sec  )  MAX  PRINCIPAL  PLANE -PHI 


I 


TIME  (sec.) 

Figure  29  -  VARIATION  OF  RMS  AND  4>L  WITH  TIME  FOR  RECORDED  ACCELERATIONS, 
STATE  BUILDING,  SAN  FRANCISCO 
CO. 5  -  SEC.  WINDOW) 
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RMS  ACCELERATION  (mm/sec  )  MAX  PRINCIPAL  PLANE  -  PHI 


RMS  ACCELERATION  (mm/sec  )  MAX  PRINCIPAL  PLANE-PHI 


In  particular,  it  was  assumed  that  the  difference  (t,)  -  t,  was  the  same 

la  1 

for  the  three  downtown  stations  in  Fig.  27.  These  estimated  values  of  t^ 
are  included  in  Table  7  and  have  been  superimposed  on  Figs.  26  through  31. 

In  general,  these  figures  show  a  behavior  similar  to  what  was  discussed  for 
the  San  Fernando  records  in  Sections  2  and  3.1:  at  t  *  t^,  the  Husid  plot 
becomes  steeper  and  there  is  an  upward  jump  in  the  RMS(t)  and  $^(t)  graphs. 
The  only  exceptions  to  this  are  the  graph  at  State  Bldg.  (Fig.  29)  and 
both  RMS  and  graphs  at  Oakland  City  Hall  (Fig.  32).  This  is  attributed 
to  the  short  time  elapsed  between  the  triggering  of  the  instruments  and  the 
first  S-wave  arrival  (0.30  sec.)  in  both  cases,  and  to  the  fact  that  a  time 
window  of  comparable  width  (6  =  0.5  sec.)  was  used  for  the  calculations. 

The  duration  of  the  strong  part  was  estimated  to  be  A  =  2  sec.  for  all 
stations,  and  the  corresponding  end  times,  t^  =  t^  +  2  have  also  been  super¬ 
imposed  on  Figs.  26  through  31. 

Figure  33  presents  a  comparison  of  acceleration  response  spectra,  com¬ 
puted  from  recorded  horizontal  accelerations  at  the  rock  site  and  at  the 
Southern  Pacific  Building  soft  soil  site.  Especially  interesting  is  the  peak 
in  the  soil  site  spectrum  at  a  period  T  =  1.1  second,  which  is  not  present 
in  the  rock  site  spectrum.  This  peak  is  due  to  the  long  period  motions  after 
t2  in  the  soil  record,  which  has  already  been  discussed  (see  Dobry  et.  al., 
1978a) . 

4 . 2  Site  Response  Simulations 

Site  response  calculations  were  performed  at  the  three  San  Francisco 
downtown  sites,  using  as  rock  input  the  horizontal  motions  recorded  at  Golden 
Gate  Park.  Before  use,  the  recorded  rock  accelerations  were  scaled  down  to 
incorporate  distance  attenuation  effects,  by  multiplying  them  by  the  same 
reduction  factors  used  by  Idriss  and  Seed  (1968),  as  follows: 
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Figure  33  -  COMPARISON  OF  RESPONSE  SPECTRA  FOR  ACCELERATIONS  RECORDED 
AT  ROCK  AND  SOIL  SITES,  SAN  FRANCISCO 
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Site 


Attenuation  Factor 


State  Building  0.80 

Alexander  Building  0.65 

Southern  Pacific  Building  0.65 

The  soil  profiles  and  geotechnical  parameters  used  for  the  simulations 
were  obtained  from  data  by  Idriss  and  Seed  (1968),  Idriss  (1969)  and  Singh 
et.  al.  (1979).  The  State  Building  soil  profile  consists  essentially  of  sand 
extending  to  the  shale  bedrock  at  a  depth  of  200  ft.,  with  occasional  layers 
of  clayey  sand.  The  Alexander  Building  soil  profile  consists  of  about  100  ft. 
of  clayey  and  silty  sand  followed  by  40  ft.  of  sand  overlying  bedrock.  The 
soil  conditions  at  the  Southern  Pacific  Building  site  are  shown  in  Fig.  34. 

The  program  SHAKE  was  used  for  the  simulations  (see  discussion  of  SHAKE 
in  Section  3.2).  In  all  calculations,  the  input  rock  motions  were  assumed  to 
act  at  a  rock  outcrop,  with  Vg  =  2,500  fps  used  for  the  San  Francisco  shale 
rock  underlying  the  three  sites.  Recent  experience  with  modulus  reduction 
curves  of  granular  soils  and  many  clays  indicates  that  the  use  of  the  stan¬ 
dard  curve  for  sand  proposed  by  Seed  and  Idriss  (1970)  may  be  generally  ap¬ 
propriate.  Therefore,  the  modulus  reduction  and  damping  curves  versus  strain 
for  sands  provided  in  SHAKE  were  used  for  the  calculations. 

Table  8  presents  a  comparison  between  computed  and  simulated  peak 

accelerations  and  their  times  of  occurrence,  (t  -  t,). 

P  1 

Figure  35  presents  comparisons  between  simulated  and  recorded  accelera¬ 
tion  response  spectra  at  the  soft  Southern  Pacific  Building  site.  Especially 
interesting  in  Fig.  35  is  the  comparison  at  a  period  of  about  1.1  seconds, 
where  the  recorded  soil  spectrum  has  a  prominent  peak.  The  peak  in  the 
simulated  curve  at  this  period  is  much  less  prominent  or  does  not  exist. 

This  contradiction  between  recorded  and  simulated  motions  computed  at  this 


66 


Depth, 

Feet 

0 


20 

55 


S  =  Unconfined 
u  Compression 
Strength 


105 

125 

135 


225 


285 


Figure  34  -  SOIL  PROFILE  AT  SOUTHERN  PACIFIC  BUILDING  SITE  USED 
FOR  SITE  RESPONSE  CALCULATIONS 
(From  Singh  et.  al.,  1979) 
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site  using  equivalent  linear  techniques  has  also  been  noted  by  Idriss  and 
Seed  (1968)  and  by  Singh  et.  al.  (1979). 

Time  domain  analyses  of  the  simulated  soil  accelerograms  at  the  three 
sites  were  performed.  Normalized  Husid  plots  and  RMS(t)  graphs  were  developed 
and  compared  with  the  corresponding  curves  recorded  at  the  sites.  The  re¬ 
sults  of  these  analyses  and  comparisons  are  presented  in  Figs.  36  through  42. 

Figures  36,  37  and  38  show  the  simulated  Husid  plots.  In  all  three 
figures,  the  Husid  plot  for  the  rock  input  (Golden  Gate  Park  recorded  motion) 
has  also  been  superimposed.  An  examination  of  the  figures  clearly  shows  that 
the  rate  of  energy  buildup  of  the  simulated  motions  is  in  all  cases  identical 
to  that  of  the  input  rock  motions  and  does  hot  predict  the  increased  durations 
and  additional  motions  occurring  after  C2  at  the  soil  sites.  These  conclu¬ 
sions  are  similar  to  those  discussed  in  Section  3.2  for  the  site  response 
simulation  at  the  Athenaeum  site  during  the  San  Fernando  earthquake.  Also, 
the  larger  contribution  to  the  Husid  plots  of  the  aftershock  acceleration  re¬ 
corded  at  the  soil  sites,  as  compared  to  rock,  is  not  well  predicted  by  the 
simulations,  as  shown  by  Figs.  36  through  38. 

Figure  39  includes  a  comparison  between  the  RMS(t)  graphs  simulated  at 
the  Southern  Pacific  soil  site  and  that  of  the  rock  input.  Two  different 
ordinate  scales  are  used  for  ease  of  comparison.  Note  that,  except  for  the 
RMS  scale  difference,  the  two  curves  in  Fig.  39  are  practically  identical, 
therefore  confirming  that  the  SHAKE  simulation  essentially  preserved  without 
modification  the  time  domain  and  duration  characteristics  of  the  rock  input 
motion. 

Figures  40  through  42  include  comparisons  of  simulated  and  recorded  RMS(t) 
graphs  at  the  three  soil  sites.  Again,  the  agreement  is  rather  poor,  both 
in  absolute  RMS  values  and  in  the  shapes  of  the  curves.  The  simulated  RMS 


70 


(t,)a=  28.44  sec. 


» 


l 


A  ^  jSU9|U  |  SDUV  |D*OJ.  4U93J3d 


CO 


0) 

l-l 

3 

00 

•H 

U-. 


71 


at  Southern  Pacific  Building  Site) 


=  0.30  sec. 


r 


1 


AIISN3J.NI  SVIHV  1V101  iN3Dy3d 


X 

< 

H 

u 

z 

o 

T3 

N 

C 

cO 

C* 

o 

i— 

X 

4J 

X 

w 

<D 

H 

r—i 

< 

cO 

X 

O 

2 

U) 

<D 

CO 

e 

•H 

Q 

4J 

z 

< 

CC  CO 

o 


K 

M 

CO 

Q 

Vh 

X 

GO 

o 

M 

O 

X 

J-i 

CO 

X 

<D 

H- 

t— s 

O 

w 

<U 

X 

H 

O 

PH 

< 

a 

H 

(0 

Q 

CO 

h-l 

6 

CO 

*■ 

O 

X 

CO 

X 

Z 

Uh 

O 

Q 

»— < 

a> 

w 

H 

u 

N 

CO 

h-i 

p2 

X 

w 

w 

X 

d) 

w 

3 

2 

u 

r—i 

o 

CJ 

CO 

z 

< 

> 

I 


m 

<U 

u 

3 

GO 

•H 

U< 


72 


recorded  at  State  Building  Site) 
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Figure  38  -  NORMALIZED  HUSID  PLOTS  OF  RECORDED  AND  SIMULATED  HORIZONTAL 
ACCELERATIONS,  ALEXANDER  BUILDING  SITE  (time  scale,  t  and 
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Figure  39  -  COMPARISON  OF  RMS  CURVES  FOR  INPUT  AND  OUTPUT  ACCELERATIONS 
SOIL  RESPONSE  SIMULATION,  SOUTHERN  PACIFIC  BUILDING  SITE 
(NOTICE  DIFFERENCE  OF  SCALES  FOR  TWO  CURVES) 

(time  scale,  t,  and  t.  values  are  from  rock  input) 


ordinates  in  the  strong  part  are  in  all  cases  significantly  larger  than  the 
corresponding  recorded  RMS  values.  The  difference  in  the  RMS  values  could  be 
reduced  if  smaller  distance  attenuation  reduction  factors  for  the  rock  input 
accelerations  were  used.  However,  the  use  of  different  reduction  factors 
most  probably  would  not  change  the  shape  of  the  simulated  wave. 

The  results  for  the  Southern  Pacific  Building  site  in  Fig.  40  indicate 
that  the  relatively  high  intensity  of  shaking  for  t  >  t^  is  not  predicted  by 
the  simulation.  Also,  the  intensity  of  shaking  of  the  aftershock  recorded 
at  the  site  is  underpredicted.  Figures  41  and  42  reveal  that  the  aftershock 
motions  predicted  by  the  simulation  at  the  other  two  downtown  sites  were  also 
lower  than  the  recorded  motions.  Idriss  and  Seed  (1968)  and  Singh  et.  al. 

(1979)  also  found  similar  results  when  using  the  equivalent  linear  technique 
for  site  response  simulations  of  the  1957  motions  at  the  Southern  Pacific 
Site.  Singh  et.  al.  have  attributed  this  effect  to  the  strongly  nonlinear 
stress-strain  behavior  of  the  soil,  which  caused  the  effective  damping  ratio 
of  the  soil  to  decrease  and  the  amplification  to  increase  after  t^. 

A  crude  approximation  of  this  nonlinear  effect  was  attempted.  In  that, 
an  additional  simulation  was  performed  at  the  Southern  Pacific  Site,  in  which 
each  of  the  two  horizontal  components  of  the  input  rock  motions  were  arbitrarily 
divided  into  two  sections:  from  0  to  4.65  seconds,  and  from  4.65  sec.  to 
the  end  of  the  record.  Each  one  of  these  input  rock  sections  was  run  separately 
through  SHAKE,  therefore  allowing  the  program  to  iterate  to  different  modulus 
reduction  and  damping  ratios  in  each  section.  The  computed  motions  at  the  top 
of  the  soil  profile  were  then  reassembled,  and  their  Husid  and  RMS  graphs  were 
calculated.  The  results  of  this  simulation  with  successive  SHAKE  runs  have 
been  superimposed  on  Figs.  36  and  40.  These  figures  show  that  the  simulation 
with  successive  runs  improves  somewhat  the  agreement  between  simulated  and 


75 


RMS  ACCELERATlOiLlmm/sec 


450 


!  f,  =  0.65  sec. 


5  10  15  20  25  30  35 

TIME  (sec.) 


Figure  40  -  COMPARISON  OF  RMS  CURVES  FOR  RECORDED  AND  SIMULATED 

HORIZONTAL  ACCELERATIONS,  SOUTHERN  PACIFIC  BUILDING  SITE 
(time  scale,  and  t^  values  are  from  recorded  accelerograms) 
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Figure  42  -  COMPARISON  OF  RMS  CURVES  FOR  RECORDED  AND  SIMULATED 
HORIZONTAL  ACCELERATIONS,  ALEXANDER  BUILDING  SITE 
(NOTICE  DIFFERENCE  OF  SCALES  FOR  TOO  CURVES) 

(time  scale,  and  t ^  values  are  from  recorded 

accelerograms) 
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recorded  aftershock  motions ,  but  does  not  change  the  poor  agreement  between 
the  predicted  and  recorded  motions  of  the  main  event  for  t  >  tj. 
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5.  SUMMARY  AND  CONCLUSIONS 

5 . 1  Summary 

The  study  presents  analyses  of  rock  and  soil  accelerations  recorded  in 
the  Pasadena  and  the  San  Francisco  Bay  Areas,  during  two  earthquakes.  These 
analyses  were  performed  using  recently  developed  time-domain  processing  tech¬ 
niques.  A  basic  objective  is  to  uncover  the  main  time  domain  characteristics 
of  strong  motion  accelerograms,  including  similarities  and  differences  existing 
between  rock  and  soil  sites. 

Four  ground  accelerograms  obtained  in  the  Pasadena  area  during  the  1971 
San  Fernando  earthquake  were  studied,  including  one  rock  site  (Caltech  Seismo- 
logical  Laboratory)  and  three  soil  sites.  Also  studied  were  five  ground 
accelerograms  obtained  in  the  Bay  Area  during  the  1957  San  Francisco  earth¬ 
quake,  including  one  station  on  rock  (Golden  Gate  Park).  The  subsurface  con¬ 
ditions  at  the  soil  recording  stations  varied  between  relatively  stiff  soils 
in  Pasadena  to  a  soft  clay  site  in  downtown  San  Francisco  (Southern  Pacific 
Building  site).  No  deep  soil  sites  in  the  sense  defined  by  Seed  et.  al.  (1974) 
were  included  in  the  study. 

The  results  of  these  analyses  of  recorded  motions  are  presented  in  Sections 

3.1  and  4.1.  The  prc-essing  techniques  used  consisted  basically  of  plotting 
graphs  of  the  variatiut  with  time  of:  a)  the  buildup  of  energy  of  the  vector 
horizontal  ground  a<releration  (normalized  Husid  plot),  b)  the  Root-Mean-Square 
vector  horizontal  ground  acceleration,  RMS(t)  for  a  moving  time  window,  6  = 

0.5  sec.,  and  c)  the  angle  between  the  principal  direction  of  the  vector  ground 
acceleration  and  the  vertical  axis,  }^(t)  for  the  same  window,  6  =  0.5  sec. 

These  graphs  are  independent  of  the  particular  orientation  of  the  horizontal 
components  of  the  recording  instruments,  and  do  not  require  any  prior  knowledge 
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of  the  fault  or  epicenter  location  with  respect  to  the  station. 

The  use  of  these  graphs  in  Sections  3.1  and  4.1  show  that  they  supple¬ 
ment  each  other  in  providing  a  better  understanding  of  the  features  of  the 
motion.  This  includes  the  determination  of  seismic  wave  arrivals  at  each 
site  and  the  consistency  of  these  arrivals  between  different  stations.  The 
graphs  are  especially  useful  in  determining  the  time  of  the  first  S-wave 
arrival,  t^.  At  t  =  t^,  the  Husid  plot  becomes  suddenly  steeper,  and  the 
RMS(t)  and  $^(t)  graphs  jump  from  low  to  high  values.  This  is  illustrated 
in  Figs.  10  and  15,  by  the  graphs  corresponding  to  the  Millikan  Library  ac¬ 
celerogram  recorded  in  Pasadena.  (When  the  jumps  do  not  occur,  as  it  is  the 
case  for  the  Oakland  City  Hall  station  illustrated  by  the  graphs  in  Fig.  32, 
this  can  be  taken  as  a  good  indication  that  the  instrument  was  triggered, 
either  after  t ^  by  the  S-wave  arrivals,  or  very  shortly  (less  than  0.5  second) 
before  t^.)  An  outstanding  example  of  this  consistent  behavior  of  the  graphs 
at  t  =  tj  is  provided  by  the  aftershock  which  occurred  in  San  Francisco  about 
30  sec.  after  the  main  event.  The  arrival  of  the  aftershock  P-waves 
is  marked  by  a  dip  in  the  <$^(t)  values  to  $  f  20°,  while  the  first  S-wave 
arrival  coincides  with  a  jump  back  to  =  80°  to  90°.  Th^s  is  true  for  the 
four  rock  and  soil  stations  which  recorded  the  aftershock,  as  shown  in  Figs. 

28  through  31. 

The  strong  part  of  earthquake  motions  at  rock  sites  has  been  defined  in 
other  publications  as  having  a  duration  A  =  t^  -  t^,  with  t^  and  t^  corres¬ 
ponding  approximately  with  the  first  and  last  S-wave  arrivals  travelling  a 
direct  path  between  the  source  and  the  site.  The  time  t^  is  correlated  with 
a  flattening  of  the  Husid  plot  and  a  decrease  in  values  of  the  RMS(t)  graph 
for  records  on  rock.  This  behavior  at  t  ■  tj  is  illustrated  for  the  two  rock 
sites  in  Figs.  9  and  12  for  the  Seismological  Laboratory  station  and  in  Figs. 
26  and  28  for  the  Golden  Gate  Park  station.  The  values  A  =  6  sec.  and 
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A  =  2  sec.,  respectively,  used  to  determine  t2  are  consistent  with  the  inde- 
pendentally  estimated  durations  of  the  rupture  at  the  generating  faults  during 
the  two  events. 

Values  of  t^  and  of  -  t^  +  A  were  superimposed  on  all  graphs  of  re¬ 
corded  motions,  thereby  cutting  the  record  in  three  parts:  i)  the  first  part 
(t  <  t^),  ii)  the  strong  part  (t^  <  t  <  t^)  and  iii)  a  third  part  (t  >  t^). 

An  examination  of  <t>^(t)  graphs  during  the  strong  part  shows  that  the  acceler¬ 
ations  were  predominantly  horizontal.  This  suggests  that  the  direct  S-wave 
arrivals  consisted  mainly  of  vertically  (or  almost  vertically)  propagating 
S-waves  (SH-waves) ,  at  both  rock  and  soil  sites  studied. 

The  separation  of  all  records  into  these  three  parts  was  especially  use¬ 
ful  for  determining  the  influence  of  site  conditions.  RMS(t)  graphs  and  Husid 
plots  of  rock  and  soil  records  were  compared  after  shifting  the  time  scales 
so  that  t^  and  t^  would  coincide  between  stations.  In  the  Pasadena  area, 
these  comparisons  show  that:  a)  a  main  wave  arrival  occurred  at  all  stations 
at  t  =  t^  +  3  (see  Fig.  18  and  t  -  t^  values  in  Table  4),  b)  the  sequence  of 
wave  arrivals  was  almost  identical  between  one  of  the  soil  sites  (Jet  Propul¬ 
sion  Laboratory)  and  the  rock  site,  without  revealing  any  significant  differ¬ 
ence  due  to  site  conditions  (see  Fig.  14),  and  c)  some  differences  between  the 
wave  arrivals  at  another  of  the  soil  sites  (Athenaeum)  and  the  rock  site  are 
observed,  including  possibly  some  additional  wave  reflections  in  the  soil  at 
t  >  t^  (see  Fig.  18).  In  the  downtown  San  Francisco  area,  where  three  soil 
stations  were  located,  significant  ground  motions  were  recorded  after  t2, 
which  were  not  present  in  the  rock  accelerogram.  These  motions,  presumably 
also  due  to  additional  wave  reflections  in  the  soil,  are  especially  prominent 
at  the  soft  S  ithern  Pacific  Building  site.  These  motions  at  the  soft  site 
have  a  period  of  about  1  second  (see  Figs.  27,  31  and  33). 


The  suggested  conclusion  discussed  above  that  the  wave  arrivals  during 
the  strong  part  of  the  motion  are  predominantly  SH-waves  is  consistent  with 
the  basic  assumption  made  when  performing  one-dimensional  site  response 
analyses.  Therefore,  site  response  simulations  were  performed  at  one  Pasa¬ 
dena  soil  site  (Athenaeum)  and  at  the  three  San  Francisco  downtown  sites. 

In  all  cases,  the  recorded  rock  motions  were  used  as  input  to  the  equivalent 
linear  code  SHAKE.  The  calculated  soil  accelerograms  were  then  processed 
using  the  time  domain  techniques  described  above,  and  the  calculated  graphs 
were  compared  with  the  corresponding  graphs  obtained  from:  a)  the  recorded 
soil  motions  at  the  same  site,  and  b)  the  rock  motions  used  as  input.  These 
results  and  comparisons  are  shown  in  Sections  3.2  and  A. 2. 

In  all  cases  studied,  the  time  domain  characteristics  of  the  simulated 
soil  motions  are  very  similar  to  those  of  the  input  rock  motions.  The 
normalized  Husid  plots  and  RMS(t)  graphs  are  about  identical  in  shape  (Figs. 
10,  23  and  36  through  39).  The  increased  motions  which  are  present  in  some 
soil  records  after  t?  are  not  well  simulated. 

A  possible  explanation  is  that  this  is  a  defficiency  of  the  equivalent 
linear  approximation  to  nonlinear  soil  behavior,  rather  than  a  basic  flaw  of 
the  one-dimensional  method  itself.  The  equivalent  linear  method  used  by  SHAKE 
defines  modulus  reduction  factors  and  damping  ratios  for  the  soil  which  are 
consistent  with  the  level  of  shaking  in  the  strong  part.  The  use  of  these 
damping  ratios  for  t  >  t2»  where  the  input  rock  accelerations  are  weaker,  may 
produce  an  underestimation  of  the  soil  motion  after  tj.  Therefore,  a  special 
calculation  was  performed  for  the  Southern  Pacific  Building  site,  using  two 
independent  and  successive  SHAKE  runs.  It  was  expected  that  this  would  pro¬ 
vide,  at  least  in  a  crude  way,  a  better  representation  of  the  nonlinear  soil 
behavior  after  t2*  However,  the  results,  shown  in  Fig.  40,  do  not  show  an  im¬ 
proved  agreement  between  simulated  and  recorded  soil  motion  after  t2* 
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5.2  Conclusions 


The  two  main  objectives  of  this  report  were:  a)  to  verify  the  use  of 
the  proposed  time  domain  processing  techniques  for  studying  recorded  rock 
and  soil  accelerograms,  and  b)  to  use  the  same  techniques  for  evaluating  the 
simulation  of  soil  accelerograms  provided  by  site  response  analyses. 

Both  objectives  were  accomplished.  Although  the  proposed  techniques 
still  require  some  further  development,  their  use  already  provides  a  simple 
and  attractive  way  of  analyzing  and  comparing  strong  motion  accelerograms. 

Specific  conclusions  for  the  recorded  and  simulated  motions  studied  are: 

1)  Husid  plots,  RMS(t)  and  ^  ( t )  graphs  show  clearly  the  beginning  of  the 
strong  part  (t  =  t^)  associated  with  the  first  S-wave  arrival,  in  rock 
and  soil  sites. 

2)  Husid  plots  and  RMS(t)  graphs  can  be  used  for  estimating  the  end  of  the 
strong  part  (t  =  t2)  of  rock  records. 

3)  The  comparison  of  normalized  Husid  plots,  and  RMS(t)  graphs  between  dif¬ 
ferent  stations,  after  shifting  their  time  axes  so  as  to  make  t^  coin¬ 
cide,  is  very  useful.  The  RMS(t)  comparisons  are  especially  fruitful, 
as  they  reveal  the  consistency  in  detail  between  wave  arrivals  at  dif¬ 
ferent  stations  during  the  strong  part,  and  the  additional  motions  pre¬ 
sent  at  some  soil  sites,  after  the  end  of  the  strong  part.  In  the  Pasa¬ 
dena  area  during  the  1971  San  Fernando  earthquake,  a  main  arrival  occurred 
in  all  stations  3  seconds  after  the  first  S-wave  arrival.  This  is  shown 
by  prominent  peaks  at  all  RMS(t)  graphs  at  t  =  t^  +  3  and  also  by  the 
times  of  occurrence  of  the  peak  accelerations  at  all  stations. 

4)  The  site  response  simulations  using  the  equivalent  linear  technique  show 
that  this  technique  computes  soil  accelerograms  which  preserve  the  shapes 
of  Husid  plots  and  RMS(t)  graphs  of  the  input  rock  motions.  This  is 
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reasonably  consistent  for  the  strong  part,  where  the  RMS(t)  graphs  of 
recorded  rock  and  soil  motions  are  generally  similar.  These  site  re¬ 
sponse  simulations  do  not,  however,  predict  the  additional  soil  motions 
observed  after  t2*  It  is  suggested  that  further  investigations,  in¬ 
cluding  truly  nonlinear  site  response  simulations  at  some  of  the  soil 
sites,  are  required  for  a  better  understanding  of  this  effect.  How¬ 
ever,  it  must  be  noted  that  nonlinear  site  response  simulations  will 
not  model  other  effects,  which  may  be  responsible  for  the  motions  after 
t2  in  some  soil  sites,  such  as  surface  waves,  soil-structures  inter¬ 
action,  etc. 
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APPENDIX  1 

VARIATION  OF  ANGLE  $  (t)  WITH  TIME 
(Units  of  are  degrees) 

RECORDED  GROUND  ACCELEROGRAMS 

SAN  FRANCISCO,  1957 
SAN  FERNANDO,  1971 
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APPENDIX  2 


VARIATION  OF  PREDOMINANT  FREQUENCY,  f(t) 
WITH  TIME 

(Units  of  f(t)  *  Time  Window  Cycles 
are  cps). 
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